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Abstract 

We discuss in this paper the phase diagram of a 
diblock polyampholyte solution in the limit of high 
ionic strength as a function of concentration and 
charge asymmetry. This system is shown to be very 
similar to solutions of so-called charged-neutral di- 
block copolymers: at zero charge asymmetry the 
solution phase separates into a polyelectrolyte com- 
plex and almost pure solvent. Above a charge 
asymmetry threshold, the copolymers are soluble 
as finite size aggregates. Scaling laws of the aggre- 
gates radius as a function of pH of the solution are 
in qualitative agreement with experiments. 

1 Introduction 

Mixtures of oppositely charged polyelectrolytes 
have been studied quite extensively over the last 
thirty years by various experimental groups jlj 0, ^| . 
There are numerous motivations for these stud- 
ies such as the similarities with biological systems 
or the applications to cosmetics and food indus- 
tries. Depending on the stoichiometry of the mix- 
ture, (the relative concentrations, the relative chain 
lengths and charge densities), one observes mainly 
two types of behavior: a macroscopic phase separa- 
tion between the solvent and the polymers, or a par- 
tial aggregation of the polymer chains || . In both 
cases, one speaks about polyelectrolyte complexa- 
tion. The electrostatic attraction between chains 
of opposite signs is clearly responsible of this be- 
havior. Despite this simple statement, there is no 
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deep understanding of the underlying mechanisms 
of polyelectrolyte complexation. The introduction 
of the polyelectrolyte multilayers concept almost 
ten years ago has revived the activity of this re- 
search area, mainly due to the numerous potential 
applications of the multilayers Q. The principle 
of the polyelectrolyte multilayers growth is to ad- 
sorb sequentially polyelectrolytes of opposite signs 
on a charged substrate. The adhesion between two 
consecutive layers is attributed to the complexation 
between the oppositely charged polyelectrolytes at 
the interface between layers where they strongly in- 
terpenetrate n ||. 

In a bulk solution, the complex formation is very 
similar to the demixing phase transition of a neutral 
polymer solution under poor solvent conditions if 
the polyelectrolyte mixture is symmetric in charge, 
i.e. if the polycations and poyanions carry the same 
charge in absolute value, and their total concentra- 
tions are the same: there is a polymer-solvent phase 
separation. If the mixture is no longer symmetric in 
charge, soluble complexes carrying a net charge can 
exist in solution H. The aggregation process de- 
pends on the charge asymmetry between the poly- 
electrolytes, but also on the relative concentrations 
of the polyions. Due to the large number of degrees 
of freedom, a comprehensive study of an asymmet- 
ric polyelectrolyte mixture is quite difficult. Ever- 
aers et al. have performed recently a study of so- 
lutions of random polyampholytes which are very 
similar systems [frl. They have allowed both for a 
macroscopic phase separation and for a finite size 
aggregation process but they have restricted them- 
selves to some specific limits for a sake of simplic- 
ity. In the case of a mixture of oppositely charged 
polyelectrolytes, one can suppress one degree of 
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freedom by considering diblock polyampholyte so- 
lutions. These copolymers are made of two oppo- 
sitely charged polyelectrolytes chemically linked by 
one end. The concentrations of polycations and 
polyanions in the solution are then equal. In par- 
ticular, each aggregate is characterized by a single 
aggregation number. In the equivalent polyelec- 
trolyte mixture, where the polyions are not bound, 
the aggregates are characterized by two aggregation 
numbers, the number of polycations and polyanions 
respectively. 

Our aim in this article is to study the phase dia- 
gram of diblock polyampholyte solutions as a func- 
tion of the copolymers charge asymmetry. This is 
a first step towards a more precise understanding 
of aggregation processes in mixtures of oppositely 
charged polyelectrolytes. Diblock polyampholyte 
solutions have been far less studied experimentally 
than the equivalent mixture of free polyions. Nev- 
ertheless, recent experiments have shown that the 
aggregates in solution are roughly spherical; they 
adsorb on charged substrate mostly as micelles, and 
not as isolated molecules H ||, [lO]. The aggregate 
sizes can be measured by Dynamic Light Scattering 
experiments and are of the same order as the size of 
an adsorbed aggregate. To our knowledge, in all of 
the existing experiments, the charged groups on the 
copolymers are pH-dependent. In particular, the 
copolymer has an isoelectric point at a value pFL of 
the pH, where the net charge vanishes. Around the 
isoelectric point, the solution is not stable: there is 
a macroscopic phase separation which is attributed 
to the complexation of the copolymers. For pH 
values not too close to pHj, aggregates of diblock 
polyampholytes form in the solution, but the solu- 
tion is macroscopically stable. Therefore, our theo- 
retical description must consider both possibilities 
of macroscopic phase separation and of partial ag- 
gregation of the chains in the supernatant. This is 
achieved in several steps. 

In a first step (section ||) Q , we neglect the influ- 
ence of the supernatant content on the properties 
of the dense phase for diblock polyampholytes with 
a quenched homogeneous charge distribution. It 
has been realized quite early that the properties of 
polyclectrolytc complexes can be described in the 
low charging limit by an equilibrium between the 
attraction induced by charge fluctuations and the 
hard core interactions : due to the very low mix- 
ing entropy of polymers, polymer mixtures are very 



sensitive to small interactions fill]. This low mix- 
ing entropy is, for example, responsible for the so- 
called "incompatibility" between neutral polymers. 
In polyelectrolyte mixtures, the phase diagram is 
determined by a subtle balance between the incom- 
patibility and complexation transitions [fl2f . Using 
our previous results on the complexation between 
oppositely charged polyelectrolytes in a symmetric 
mixture, we describe the swelling of the dense phase 
by increase of the counterion excess osmotic pres- 
sure with the asymmetry of the copolymers 0. In a 
second step (cf section ||), the aggregation process is 
described in the absence of any macroscopic phase 
separation. By analogy with the associating be- 
havior of charged-neutral diblock copolymers (ion- 
ized hydrophilic block-neutral hydrophobic block), 
we propose that the copolymers form spherical mi- 
celles with a neutral core made of symmetric poly- 
electrolyte complex, and branches carrying the net 
charge of the aggregates. Finally, the phase di- 
agram is drawn by combining all these results in 
section |I| We also discuss in this section the ef- 
fect of "charge annealing" on the properties of the 
aggregates in the case where the ionizable groups 
along the blocks are weak acids and weak bases. 
All the energies in this paper are measured in units 
of the thermal excitation ksT. 

2 Complexation thermody- 
namics of an asymmetric 
polyelectrolyte mixture 

We consider in this section a monodisperse solu- 
tion of diblock polyampholytes of total degree of 
polymerization N. We assume that all the chem- 
ical properties of the blocks are identical except 
their charge that are of opposite signs. Let N + , f + 
and AL. , /_ be respectively the degrees of poly- 
merization and charge densities of the polycationic 
and polyanionic blocks of the copolymer. We ne- 
glect in this paper the effects of charge distribution 
and assume that the total charge of each block is 
smeared out along chain. We assume that the sol- 
vent is a O-solvent and we treat it as a dielectric 
continuum of dielectric constant e. The monomer 
size is denoted by a. It is also convenient to in- 
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troduce the Bjerrum length Is = ^tk B T > which 
represents the typical distance between two ele- 
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mentary charges having an electrostatic energy of 
ksT. When the polyampholyte is symmetric in 
charge, i.e. N + f + = N-f-, the solution phase sep- 
arates into a dense phase containing almost all the 
copolymers and a dilute phase composed mainly 
of solvent. Within the dense phase, the diblock 
copolymer concentration (expressed as a number of 
chains per unit volume) Cdense is related to cationic 
and anionic monomer concentrations c + and c_ by 
c+ + c_ = c dense N. When N+f+ ^ N-f-, there 
is an excess charge in the dense phase, which is 
neutralized by the small ions of the solution. The 
difference in translational entropy of these ions be- 
tween the two phases produces an excess osmotic 
pressure which tends to swell the dense phase. Fol- 
lowing the results of ref |J , it is possible to describe 
the effect of a charge excess in the dense phase by 
a generalized two phases model. The dense phase 
is mainly composed of polyampholytes, while the 
dilute phase is modeled as a simple electrolyte. 
This amounts to neglecting the formation of finite 
size aggregates in the supernatant of the solution. 
This approximation is discussed in the next section. 
The concentrations of monovalent small ions in the 
dense phase and in the dilute phase are respectively 
n± and n±o = uq. The conditions for equilibrium 
between the two phases are given by the equality 
of various small ions chemical potentials in the two 
phases, and the equality of the osmotic pressures. 
Therefore, we need to compute the free energies in 
the two phases to write the equilibrium conditions. 

As previously mentioned, the formation of poly- 
electrolyte complexes in symmetric polyclectrolyte 
mixtures cannot be described at the simplest mean 
field level: because of the global electroneutrality, 
there is no pure coulombic term in the mean field 
free energy. One has to include at least fluctuations 
around the electroneutral state at the lowest level 
to describe the complexation. This can be done 
within the RPA formalism [[t3| [l4| . In the case of a 
diblock polyampholyte solution, we will assume in 
a first approximation that the RPA correction to 
the mean field free energy of the dense phase is the 
same as that of the equivalent polyelectrolyte mix- 
ture where the two blocks are not linked: our aim 
in this paper is just to describe qualitatively the 
phase diagram at the level of scaling laws. Never- 
theless, it seems possible to calculate more precisely 
the RPA free energy of the copolymers by using the 
results of reference [|l2|. This will be the topic of 



a future work. The RPA free energy density of the 
dense phase is written 
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The first term of eq. ([j]) is the third virial contri- 
bution. This is the first hard core repulsive inter- 
action for a polymer in a O-solvent. This contribu- 
tion is necessary to insure the stability of the com- 
plex. The second term is the translational entropy 
of small ions in the dense phase. Notice that the 
translational entropy of the copolymers has been 
neglected because it is of order . Finally the last 
term is the RPA correction to the mean field free 
energy. It is calculated by associating a Gaussian 
statistical weight to the concentration fluctuations 
neglected in the mean field theory [|l3| . In the case 
of a symmetric mixture (/+c + = /_c_), the exact 
result reads 
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We introduced in the last equation three character- 
istic lengths £ w , q+, K _1 which are defined as fol- 
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k 2 = AttIb (J2i n i)- The first length is the cor- 
relation length of the equivalent neutral polymer 
mixture in a O-solvent. The two other lengths are 
associated respectively to the screening of electro- 
static interactions by the polyelectrolytes and by 

the small ions. The dimensionless ratio of these 
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two lengths is denoted s = t-. Notice that in the 
symmetrical case, the contribution of fluctuations 
associated respectively to the third virial term and 
the electrostatics are not coupled. The free energy 
density of the dilute phase is simply given within 
the RPA by the Debye-Huckel theory 111 
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with Kq = 4:7tIb n io)- I n the limit of high ionic 
strength, one can solve the equilibrium equations 
using an expansion in the small parameter s$ . Ad- 
ditionally, an expansion up to first order in kqI b en- 
sures that the Debye-Hiickel approach is valid. In 
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the symmetric case, the osmotic balance between trivial order 
the two phases reads the lowest non-trivial order 
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The effect of fluctuations associated to the third 
virial term can be absorbed by a renormalization 
of w. The renormalized third virial term is denoted 
w. This osmotic balance lead us to introduce the 
notion of 'complexation blobs' [pL 16 : by analogy 



with neutral polymers under poor solvent condi- 
tions, the dense phase can be viewed as a compact 
packing of complexation blobs of size £ c ~ Kq 1 Sq- 
At length scales below £ C) the concentration fluctu- 
ations induced attraction is not relevant compared 
to the thermal energy fc^T, while it is dominant 
at larger length scales. This interpretation of the 
structure of the dense phase will be used in the next 
section to discuss the formation aggregates. 

For asymmetric mixtures (/+c+ =/= /_c_), the 
calculation of AFdense is more complicated. Never- 
theless we are looking for the effect of a small charge 
asymmetry on the structure of the dense phase. 
The calculation can be therefore simplified by ex- 
panding the RPA correction term with the charge 
asymmetry [(L6) . The main difference with the sym- 
metric calculation of AF ( i en se is the appearance of 
a new term coupling the fluctuations associated to 
the third virial coefficient and the electrostatic in- 
teractions. It can be shown a posteriori that this 
term simply renormalizes the prefactors in the os- 
motic balance. It will be therefore neglected. The 
equilibrium of the small ion chemical potentials be- 
tween the two phases requires the introduction of 
a Donnan potential, because of the electroneutral- 
ity constraint in the two phases [jlTj . Therefore the 
concentrations of small ions in the dense phase are 
related by 
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with the chemical potential difference due to the 
polarization energy (RPA correction term) A/i po / = 



AFdUute)- This equation in n + can 



be solved with the constraint of electroneutrality at 
the first order in kqIb- Finally, the osmotic balance 
between the two phases reads at the lowest non- 
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We used in the preceding equation the definition 
of the monomer concentrations in the dense phase 
c± = CdenseN±. The first and third terms of this 
equation are reminiscent of the osmotic pressure 
of a symmetrical complex. The second term is of 
purely entropic origin: the difference in transla- 
tional entropy of the small ions between the two 
phases induces an excess osmotic pressure inside 
the dense phase. The excess charge per copolymer 
is conveniently measured by the introduction of the 
average net charge per monomers 
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This is the net charge of a copolymer redistributed 
on all the monomers. In the limit of low asym- 
metry, the third virial term dominates over the os- 
motic term. The copolymer concentration inside 
the dense phase is therefore given by its value in 
the symmetrical case 
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The number of monomers inside a complexation 
blob reads 
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The size of the complexation blob is given by £ c ~ 
<7c a(w^/ s /a). Neglecting all numerical prefactors, 
a careful analysis of eq. @ shows that the third 
virial contribution dominates over the osmotic one 
if 
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This inequality is valid as long as | /+ — /_ 1 1 ^ L a — < 
(Z+J-) 1 / 2 . It means that the properties of the 
dense phase are well described by equations (^) and 
S 1, where / is the geometric average of 
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charge densities. If J? > 1, the osmotic term domi- 
nates over the third virial term, and the copolymer 
concentration of the dense phase can be calculated 
using the osmotic balance eq. (||). However, as will 
be shown in the next sections, finite size aggregates 
are more favorable in this range of parameters. 



3 Micelle formation 

We propose in this section a model for the aggre- 
gates in the supernatant. It has been shown in 
the last section that the polymeric net charge in 
the dense phase induces an excess osmotic pres- 
sure which tends to swell this phase. By analogy 
with the associating behavior of diblock copolymers 
charged hydrophilic/hydrophobic, we will assume 
that spherical micelles are formed in solution []l8|| : 
the polymeric net charge of the core vanishes, and 
the excess charge is localized on the surrounding 
arms in the corona of the micelle (cf figure [l]) . For 




Figure 1: Copolymer micelle. The circles of grow- 
ing diameters are the compact packing blobs of the 
Daoud Cotton model. For the sake of simplicity, 
only charges of the arms are explicitly sketched. 



a positive excess charge per copolymer, the num- 
ber of cationic monomers in the core N c+ is given 
by the electroneutrality condition N c+ f + = 7V_/_. 
The total number of monomers per chain in the core 
is denoted by N c — N c+ + N- . The excess charge 
of the corona per copolymer is Nq. In this ge- 
ometry, the non-compensated charges of the arms 



are further away from each other than in a sin- 
gle homogeneous dense phase. It is also the case 
for charged-neutral diblock copolymers where this 
geometry has been observed experimentally JOl. 
When the charge asymmetry is small enough, we 
expect other aggregate geometries to be more fa- 
vorable than the spherical one, such as cylindrical 
micelles or lamellae (planar geometry) p(|. How- 
ever it can be shown a posteriori that those geome- 
tries are less favorable at the level of scaling laws 
than a macroscopic phase separation for low charge 
asymmetry fl6|| . Therefore we will not consider ex- 
plicitly those geometries further. Before drawing 
the full phase diagram of the solution, we first con- 
sider in this section, the properties of the optimal 
spherical micelle, which minimizes the free energy 
per copolymer. We suppose therefore that all the 
required conditions for the formation of spherical 
micelles are met. The main contributions to the 
micelle free energy per copolymer are the surface 
tension of the core and the electrostatic free energy 
of the corona. In the range where spherical micelles 
are favorable, the stretching energy of the chains in 
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The "neutral" core of the micelle can be viewed 
as a symmetrical polyelectrolyte complex. Using 
the interpretation of its structure in terms of com- 
plexation blobs in the limit of high ionic strength, 
its negative free energy is evaluated with the ansatz 
l ksT per blob'. Similarly, the interfacial tension 
between the complex and the surrounding solvent 
is estimated by associating an energy fc^T per blob 
on the surface of the core. Therefore the surface 
tension reads 7a ~ ^ ~ For an aggre- 

gate containing p chains, the core size R c is given 
by ~ pN c gV 2 a? . The interfacial free energy per 
copolymer reads 
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In the limit of high ionic strength n Q a 3 > 
it has been shown that the corona of a 



.4/3 
J + 



(In/a) 1 ' 3 ' 

charged micelle can be described by the Daoud- 
Cotton model for neutral star branched polymers 
p^ . This comes from the fact that in the high 
ionic strength limit, the electrostatic interactions 
are strongly screened. This effective short range 
interaction is equivalent to an excluded volume in- 
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teraction with an electrostatic excluded volume pa- 

f 2 

rameter v e i ~ ~ u e ;a 3 . According to the Daoud- 
Cotton model, the star is composed of a compact 
packing of correlation blobs of size £ (r), where r 
is the distance to the center of the star p3[ . This 
condition implies £ (r) ~ . The chain statis- 
tics of inside a blob the excluded volume statis- 
tics £ (r) ~ 5 (r) 3 ^ 5 Ug/ 5 a 2 / 5 . The monomeric con- 
centration profile of the corona is then given by 



c (r) a 3 
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The extension of the arms 



and the free energy per chain are calculated by 
counting the number of blobs, or equivalently by 
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N c+ be the number of monomers 



Let N b = N 
per branch in the corona. In the limit R/R c >> 1 
where the spherical aggregate is the most favor- 
able geometry, one finds for the corona extension 
the same result as in a Flory-like mean field eval- 
uation balancing the osmotic pressure due to the 
salt in the micelle and the stretching of the arms 
R ~ aAT 3/ V /5 v e Y 5 - The Daoud-Cotton model 
takes into account excluded volume correlations, 
and therefore leads to a more accurate result for 
the free energy than the mean-field prediction. The 
free energy per copolymer reads 
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It depends only very weakly on Nb (logarithmi- 
cally). 

The minimization of the total free energy per 
copolymer leads to the following radius, aggrega- 
tion number and free energy for the spherical mi- 
celle 
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It can be shown that the spherical geometry is fa- 
vored if one ignores macroscopic phase separation 
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if >> c rfk — • This condition is also equiv- 
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alent to R/R c » 1. 

4 Phase diagram 

Our aim in this section is to link the the results 
of the two previous sections: in section ^ we have 
studied the macroscopic phase separation of the so- 
lution by assuming a macroscopic phase separation 
ignoring the finite size aggregation in the super- 
natant; in section Q we have evaluated the opti- 
mal properties of a spherical aggregate, neglecting 
any concentration effect. In fact, a diblock polyam- 
pholyte solution should be fully described by a 
model allowing both for macroscopic and meso- 
scopic phase separations. Let <fi, Cdense, n+ and n_ 
be respectively the volume fraction, the copolymer 
concentration, the positive and negative small ion 
concentrations in the dense phase of a macroscop- 
ically separated solution. In the dilute phase, we 
define respectively c p ,n + o and n_o the concentra- 
tions of p chains aggregates, of positive and nega- 
tive small ion. The free energy density of the solu- 
tion reads 

F — (f)Fdense ((-dense : ^-f ; ) 

+ (1 - (j>) F d ii ({c p } ,7i +0 ,n_o) (18) 
with the free energy density of the dense phase 

Rdense — Cdense (log Cdense 1) 

+ Y / n l (\ogn l -l) + F D (19) 
i=± 

and the free energy density of the dilute phase 
Fan = ^2 c p (l°g c p ~ 1 + Fp) 

p=l,2... 

+ (log n i0 -1) + F d0 (20) 

i=± 

The quantities Fjj , F p and Fdo are respectively the 
free energy of the dense phase, and the free energies 
of an aggregate of p chains and the Debye Hiickel 
polarization energy of the dilute phase. Within the 
RPA, Fd includes the excluded volume contribu- 
tion and the polarization energy of the dense phase 
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cq. (0), and depends on the various concentrations. 
Finally, F p has been evaluated in section |^. We will 
keep these very general notations in the following, 
since they are not model dependent. The equilib- 
rium between the two phases is obtained by mini- 
mizing the free energy density with respect to the 
volume of the dense phase and the various concen- 
trations, with the constraints of mass conservation 
of copolymers and small ions, and electroneutrality 
of both phases. The minimization with respect to 
the volume of the dense phase leads to the osmotic 
balance. When the contribution of translational 
entropy and polarization energy of the aggregates 
are neglected, this leads to the analysis of section 
||. This approximation is justified since the transla- 
tional entropy of an aggregate of p chains is of order 
^ , and therefore the contribution to the polariza- 
tion energy is also negligible. The other equilibrium 
equations read 

n+ = n +0 exp (-*) cxp (- (p t - ^ i0 )) 

n- = n-o cxp (+^) cxp (- (m - fMo)) (21) 



= ^dense ex P (+pNqV) exp (- (F p - pud)) 



The chemical potentials /i Q in the last equations are 
obtained by derivation of the free energy densities 
with respect to the corresponding concentration; 
they do not include the contributions of the trans- 
lational entropies; the subscript o corresponds to 
the dilute phase. The electroneutrality constraint 
in the two phases induce an electrostatic potential 
difference between the two phases 1 i>: this is the 
so-called Donnan potential. By using the first two 
equations of (^l]), we can calculate ^ at first order 
in k I b 
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The last equation of (|21|) can be rewritten c p = 
Cdense exp (— f2 p ) with the thermodynamic grand 
potential 

Q p =F p ~p (JJ, D + log C dense + Nq^) (23) 

The total copolymer concentration c tot is written 
as 
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(Y^pPCdense exp (-fip)) 
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This equation shows that the contribution of ag- 
gregates is negligible compared to the contribution 
of the dense phase if fl p > 0. In the opposite case 
ftp < 0, the concentration of aggregates becomes 
relevant. Therefore we can write a simple solubil- 
ity criterion: the chains are soluble as aggregates 
when 



F p 

-=- < fJ-D 
P 



Cdt 
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(25) 



We used in this equality the average aggregation 
number p = p that minimizes the free energy per 
diblock copolymer. This inequality means that 
copolymers are soluble when the chain chemical po- 
tential in an aggregate is smaller than the chain 
chemical potential in the dense phase. The dense 
phase chemical potential has two main contribu- 
tions: the attractive energy, — fcgT per blob, and 
the energy associated to the loss of translational 
entropy of the small ions. A careful analysis of the 
inequality (^) leads to the simplified criterion 



(Nq) 2 



. Cdense 
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Therefore the chains are soluble when 



q>f- 



.,1/3 



(26) 



(27) 



(24) 



(Note that the prefactor on the right hand side 
must be smaller than one for consistency). At the 
level of the scaling laws, we found that the dense 
phase disappears when % > 1, where / is the av- 
erage charge density of the copolymer. The chains 
are thus soluble as spherical aggregates when the 
excess charge per monomer is of the same order of 
the charge density that induces the collapse of the 
complex by fluctuations. This value corresponds at 
the level of scaling laws to the regime where the 
excess osmotic pressure due to the small ions be- 
comes of the order of the third virial contribution 
in the analysis of the osmotic balance in section |[ 
The criterion ( p7| ) does not depend on concentra- 
tion effects. This means that the solution has to be 
already phase separated or partially aggregated for 
this criterion to hold. It happens only above some 
critical copolymer concentration. 

To include concentration effects in this descrip- 
tion, one has to consider the two following asymp- 
totic limits: a micellization with no macroscopic 
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phase separation, and a macroscopic phase sepa- 
ration with no micellization. This allows to com- 
pute respectively the critical micellar concentration 
(CMC) and the critical macroscopic phase separa- 
tion concentration, and to compare them. In the 
case of a pure micellization, the total concentration 
of copolymers is written as 



PCp =C 



(28) 



The equilibrium conditions between all aggregates 
can be written as 



cf exp [- (F p - pFi 



(29) 



where F\ is the free energy of one isolated copoly- 
mer in solution. This energy contains mainly two 
contributions: the interfacial tension between the 

/ \ 2/3 

complexed part and the solvent ( J ; an d the 



attractive energy 
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Up to entropic effects, mi- 



celles are formed in solution if F% > -# , where p = p 
minimizes the free energy per copolymer in an ag- 
gregate. The CMC is then evaluated by 



log c„ 



-[Fx 



El 
p 



(30) 



For spherical micelles, this concentration can be 
rewritten as 

l...v„„- I (— ) 23 -f— V) (31) 




9v cl 



This CMC is exponentially small, as for any diblock 
copolymer system. 

In the limit of a pure macroscopic phase sepa- 
ration, with no micellization, the total copolymer 
concentration reads 



<j>c dense + (1 - 4>) a = c f 



(32) 



The equality of diblock copolymer chemical poten- 
tials between the two phases leads to 



ci = c de 



,exp[-(i?i-(Mfl+JV 9 *))] (33) 



We introduced in the last equation the Donnan po- 
tential that takes into account the loss of transla- 
tional entropy of small ion in the dense phase. The 
critical concentration for macroscopic phase sepa- 
ration is thus given by 



logc 



sep 



(Fi - ((Md + NqV))) (34) 



This critical concentration is approximately of the 
same order of the CMC. More precisely, the ratio 
between those two critical concentrations reads 



log- 



-■sep 



P 



(35) 



The right hand side of the preceding equation gives 
our simplified criterion for the solubility of aggre- 
gates. Therefore, we can draw qualitatively the fol- 









q 



Figure 2: Phase diagram (c*°* versus q) in the limit 
of high ionic strength. 



lowing phase diagram (cf figure ||) : the macroscopic 
phase separation takes place first when % < 1 and 
c tot > c sep . For higher asymmetries j > 1, spher- 
ical micelles are formed when c tot > c mc . We 
thus predict that the only aggregates to be observ- 
able when the solution is macroscopically stable are 
spherical micelles. This has to be compared to the 
recent experiments on diblock polyampholyte so- 
lutions |, As was mentioned in the 
introduction, one observes indeed in those exper- 
iments spherical aggregates in solution: their ra- 
dius is measured by Dynamic Light Scattering or by 
AFM experiments, since those copolymers mainly 
adsorb on charged substrates as micelles and not as 
single molecules. 

The only available experiments on diblock 
polyampholytes to our knowledge are performed 
with weak acidic and basic groups on the two 
blocks. It means therefore that the excess charge 
q depends on the local pH of the solution p5| . In 
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particular, the radius of the aggregates varies with 
the pH. As we shall see below in an appendix, 
our model predicts qualitatively the same varia- 
tion if local pH effects are neglected. By local 
pH effects, we mean inhomogeneities of the solu- 
tion (such as aggregation) that lead to inhomo- 
geneities of the pH that can take locally a value 
different than the imposed value. For example, the 
pH inside an isolated star-branched weak polyelec- 
trolyte is different from the value in the bulk solu- 
tion 0| . Neglecting local pH effects, we assume in 
a first approximation, that pH variations in a so- 
lution of diblock polyampholytes change only the 
value of q. Notice that in the most general case, 
N c and Nf, depend in a complex manner on q. It is 
possible to simplify those expressions for two lim- 
iting cases of asymmetry (7V + = N and / + ^ /_) 
and ( N + ^ iV_ and /+ = /_ ) . We will not analyze 
those particular cases in this paper. For the sake 
of simplicity, we consider N c and Nb as constants. 
The charge densities along the two blocks are given 
by the following action mass laws 

t^t = — (36) 

T^V = — = K B 10 U c H+ (37) 

1 - /+ C H- 

with the dissociation constants Ka and Kb of 
the acidic and basic groups. The pH of the so- 
lution is related to the concentration of H + ions by 
Ch+ = lQ~ pH ■ In cq. (|37j), the ionization product 
of water has been used. For weakly charged poly- 
electrolytes (/+,/- << f), the product of charge 
densities /+/- is approximately constant and does 
not depend strongly on the pH. As it was shown in 
the first section of this paper, the properties of the 
dense phase (in macroscopically phase separated 
sample or in the core of a micelle) depend only on 
/+/- (for example g c ~ (/+/-)~ 2 )- Therefore pH 
variations do no affect strongly the core properties 
of spherical micelles within this very crude model. 
The only effect of pH variations is to change the 
corona extension which depends strongly on the 
charge density of the arms. For a positive excess 
charge (low pH values) , the radius of a micelle vary 
as 

2/5 



charged micelle (high pH values) scales like 



R 



(38) 



R ■ 



-2/5 



-H+ 



(39) 



Therefore the extension of the micelle increases 
with the pH. For intermediate values close to the 
isoelectric point pH^, the excess charge vanishes 
and the solution macroscopically phase separates. 
Notice that eq. ( |f5|) is no longer valid when R ~ R c 
since the logarithmic factor becomes important. In 
this range of parameters, the solution has already 
phase separated as it discussed at the beginning of 
this section. Close to this regime the aggregation 
number increases again with the pH for pH<pH; 
and decreases for pH>pHi. 

The qualitative behavior of the micelle radius 
with pH is qualitatively observed in experiments of 
reference ||. The authors study the adsorption of 
dilute diblock polyampholyte solutions. As it pre- 
viously mentioned, AFM experiments on the sub- 
strates and Dynamic Light scattering experiments 
in solution confirm the presence of spherical aggre- 
gates. By studying the adsorption kinetics as a 
function of the ionic strength of the solution, a dif- 
fusion coefficient towards the surface is measured 
and can be used as an alternative mean to deter- 
mine the size of the aggregates. This radius is found 
to decrease with added salt, in qualitative agree- 
ment with our model which predicts R ~ n Q 13 / 25 . 
Notice that no significant variations of the aggre- 
gate radius with pH or ionic strength is observed 
in the experiments of Goloub et al. with another 
kind of diblock polyampholytes. We believe that 
this is probably due to the relevance of non electro- 
static interactions in this case 10 . It seems thus 



The radius of a spherical aggregate decreases with 
the pH. On the contrary, the radius for a negatively 



that diblock polyampholyte aggregates can be de- 
scribed by a very crude model of polyampholytes 
with quenched charge distribution with no local pH 
effects. The agreement with the experiments is not 
quantitative, but all the qualitative trends are pre- 
dicted. We give some arguments in the appendix 
of this paper that justify partially this very crude 
model. 



5 Concluding remarks 

We have studied in this paper the phase diagram of 
monodisperse diblock polyampholyte solutions. In 
the case of charge symmetry (N + f + = N-f-), the 
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properties of the solution are very similar to those 
of an equivalent homopolymcr mixture obtained by 
cutting the junction point between the oppositely 
charged polyelcctrolytcs in the copolymer : there is 
a macroscopic complexation phase separation be- 
tween the polymers and the solvent. The structure 
of the dense phase has been here described in terms 
of complexation blobs. When the charge asym- 
metry increases, the small ions of the solution are 
strongly coupled to the copolymers: the electroneu- 
trality of the dense phase requires a finite density of 
countcrions that induces an excess osmotic pressure 
which tends to swell the complex. Above a charge 
asymmetry threshold, this excess osmotic pressure 
is larger than the attractive energy collapsing the 
complex: the copolymers are soluble in the form 
of finite size aggregates. We propose in this paper 
that the chains are aggregated in solution as spher- 
ical micelles with a neutral core, the excess charge 
being carried by the arms of the corona. In the 
limit of high ionic strength, the micelle corona is 
described by the Daoud Cotton model for neutral 
polymeric stars with an effective electrostatic ex- 
cluded volume v e i = 47r a . The structure of the 

K o 

micelle core is very similar to a neutral polymer 
solution in a poor solvent and is described using 
complexation blobs. 

The only available experiments that we are 
aware of on diblock polyampholyte solutions are 
performed with annealed polyampholytes carrying 
acidic and basic groups along each blocks. There- 
fore the charge distribution along the chain is an- 
nealed and depend on the local pH of the solution. 
By neglecting any local pH effects, it is possible to 
predict at least qualitatively the variation of the 
micelle radius with pH and salt amount. In the 
limit of high ionic strength, this simple model can 
be partially justified: the charge fluctuations still 
collapse the core of the micelles and are in fact 
stronger in the case of an annealed charge distri- 
bution. Moreover the pH variations of the buffer 
are transmitted to the corona and the core of the 
micelles. Of course this very simple model is not 
expected to be a quantitative description of the ag- 
gregates because of the naive modelization of the 
annealed effect. The main problem for a more pre- 
cise modelization comes from the dependencies of 
some parameters of the complex as a function of 
the charge asymmetry: the number of monomers 



in the core and in the corona of the micelle, the 
charge densities inside the core and on the arms of 
the corona, etc... 

The micellization model proposed in this paper 
is a first step towards a more precise understanding 
of asymmetries effects in mixtures of (unbound) op- 
positely charged polyelectrolytes. In this case, the 
separation of the polyions produces a new degree of 
freedom. In particular, the aggregates are charac- 
terized by two aggregation numbers. This will the 
topic of future work. 
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Appendix A 

We give briefly in this appendix some arguments 
justifying the very crude model of charge annealing 
in the micelles in the limit of high ionic strength. 

A.l RPA free energy 

To understand the effect of an annealed charge dis- 
tribution, let us see how the polarization energy 
driving the complexation is modified. We consider 
Gaussian chains for the sake of simplicity (no third 
virial contribution). The general form of the RPA 
correction term to the free energy is given by 



dot 



a- 1 

'-'RPA 

2irV 



(40) 



where G~^ PA is the inverse structure matrix matrix 
at the RPA level jl3) . The determinant of this ma- 
trix in the case of an annealed charge distribution 
is exactly the same as in the quenched distribution 
case if the Debye-Hiickel screening length is defined 
as 



k 2 = AttIb n % + /+ (1 - /+) c+ + /_ (1 - /_) c_ 

(41) 

The mobile charges on the chains contribute to 
the screening of electrostatic interactions. The 
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quantities rij are the concentrations of small ions 
(salt ions, H + ,OH~). Notice that we will some- 
times denote the concentrations of H + and OH~ 
by n-H+ = c H + and n OH - = c OH - . Therefore the 
RPA correction term reads 



AF 



1/2 



(42) 



48^ B (/3 C+ +/2 C _) 



and s = % . It is ex- 

li 



with qf 

actly the same as for quenched distributions, except 
for the new definition of n. The mean field contri- 
bution to the free energy density is given by 

_F_ C+ + C- 

kf 



-(log- 



N y ° N 
^2m(logrii - 1) 



1) 



+c_[/_log/_ + (l-/_)log(l-/_) 

+/-M-A- + (1- /-)/Mh] 

+c+[/+ log/+ + (1 - /+) log(l - /+) 

+ f+fl B+ + (1 - f+^BOH] (43) 

(44) 

where we take into account the translational en- 
tropies of the chains and the small ions in the so- 
lution, the mixing entropy of the small ions on the 
chains and their equilibrium with a charge reser- 
voir. 

A. 2 Equilibrium between two phases 

As in section |^, we neglect the influence of copoly- 
mers in the dilute phase. We can write the equi- 
librium conditions as in section ^|. In the limit of 
high ionic strength, one can solve at the first order 
in kqIb the equilibrium equations for the small ions 
concentrations fl6| . This leads to introduce a new 
parameter 77 = /+ (1 - /+) c+ + /_ (1 - /_) c_. In 
the symmetric case where /+c+ = /_ c_, the os- 
motic balance reads at lowest order in kqIbjSq 1 
and^- 

n 



An = 



24n 1 



3 V 



(45) 



2 2n Q u ' 

Thus we find that at lowest order that charge an- 
nealing leads to an increased attractive term free 
energy the complex: this is quite natural since there 
is a new fluctuating variable. Nevertheless, the 
dominant term is still the same as in the quenched 
case. 



A. 3 Local pH effects 

We denote in this section with a "0" subscript the 
quantities related to the region outside the micelles, 
namely the bulk solvent properties. In particular, 
the charge density of an isolated weak polyelec- 
trolyte (no local pH effects) is fa. In the case of 
star-branched weak polyelcctrolytes, it was shown 
that the charge density of the arms depends on the 
pH of the buffer solution like |22] 

f 2 = — (cr+O + n sel + n strong base) (46) 

if /o_c >> (ctf+o + n se i + n strong base ). Notice that 
the contribution of the polarization energy has been 
neglected. In the opposite case, the charge density 
is / = fo- It means therefore that in the high ionic 
strength limit, the pH variations in the solution 
buffer are entirely transmitted to the inside of the 
star. 

All the arguments presented above can still be 
used: in the limit of high ionic strength, local pH 
effects are very weak, and therefore they can be 
ignored in first step of modelization. 
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